function [cust_cp] = custom_temporal(xx,t)
    %return cust_cp with the value of Cp for each node
    [nnod,ndm] = size(xx);
    cust_cp = zeros(nnod,1);
    %% completar desde aqui
end